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Abstract. Dynamo action is investigated in simulations 
of locally isotropic and homogeneous turbulence in a slab 
between open boundaries. It is found that a "pseudo- 
vacuum" boundary condition (where the field is verti- 
cal) leads to strong helicity fluxes which significantly re- 
duce the amplitude of the resulting large scale field. On 
the other hand, if there is a conducting halo outside the 
dynamo-active region the large scale field amplitude can 
reach larger values, but the time scale after which this field 
is reached increases linearly with the magnetic Reynolds 
number. In both cases most of the helicity flux is found 
to occur on large scales. From the variety of models con- 
sidered we conclude that open boundaries tend to lower 
the saturation field strength compared to the case with 
periodic boundaries. The rate at which this lower satura- 
tion field strength is attained is roughly independent of 
the strength of the turbulence and of the boundary condi- 
tions. For dynamos with less helicity, however, significant 
field strengths could be reached in a shorter time. 

Key words: MHD - Turbulence 



1. Introduction 

Significant progress has been made in recent years in the 
understanding of the generation of large scale magnetic 
fields in astrophysical bodies by dynamo action. Some es- 
sential features of mean-field dynamo theory (e.g., Mof- 
fatt 1978, Krause & Radler 1980) are recovered qualita- 
tively and quantitatively in direct three-dimensional sim- 
ulations, but the time scale on which large scale fields are 
established becomes very long as the magnetic Reynolds 
number increases (Brandenburg 2001a, hereafter referred 
to as B2001). This result, which has been obtained us- 
ing simulations of helical turbulence in a periodic domain, 
is now understood to be a consequence of the fact that 
the large scale magnetic fields generated by such flows 
possess magnetic helicity which is conserved, except for 
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(microscopic) resistive effects. This behaviour can also 
be reproduced with a mean-field model with resistively 
dominated ('catastrophic') quenching of a-cffect and tur- 
bulent magnetic diffusivity; see B2001 for details. The 
implications of boundary conditions for 'catastrophic' a- 
quenching were first pointed out by Blackman & Field 
(2000a), who showed that the results of earlier simula- 
tions by Cattaneo & Hughes (1996) could be understood 
as a consequence of periodic boundary conditions. 

The problem is that in order to generate large scale 
fields with finite magnetic helicity on time scales shorter 
than the resistive time, one also needs to generate mag- 
netic fields with the opposite sign of magnetic helicity, 
so that the net magnetic helicity stays close to its ini- 
tial value. In an astrophysical body the sign of the mag- 
netic helicity is different in the two hemispheres, but pre- 
liminary models with helicity flux through the equatorial 
plane suggest that this actually diminishes the net large 
scale magnetic field rather than enhancing it (Branden- 
burg 2001b). The other possibility is that one can have a 
segregation in scales (e.g. Seehafer 1996) such that large 
scale magnetic fields with magnetic helicity of one sign 
are generated within the domain of interest, whilst small 
scale magnetic fields with opposite sign of magnetic he- 
licity are expelled through the boundaries (Blackman & 
Field 2000b) . Whether or not this mechanism really works 
is unclear. So far there has only been a phenomenological 
approach suggesting that the a-effect would be large if a 
flux of magnetic helicity through the boundaries was pos- 
sible (Kleeorin et al. 2000). 

The aim of the present paper is to investigate the effect 
of open boundaries on the dynamo. A simple local bound- 
ary condition is the vertical field condition (sometimes also 
referred to as pseudo- vacuum boundary condition), that 
was used previously both in magnetoconvection (Hurlburt 
& Toomre 1988) and in accretion disc simulations showing 
dynamo action (Brandenburg et al. 1995). This boundary 
condition allows for a finite magnetic helicity flux through 
the boundaries, even though the flux of magnetic energy 
(i.e. the Poynting flux) is still zero. We also consider cases 
with a conducting halo outside the disc, as well as cases 
with density stratification, so that magnetic buoyancy can 
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contribute to turbulent flux transport out of the dynamo- 
active region. 

2. The model 

We solve the isothermal compressible MHD equations for 
the logarithmic density In p, the velocity it, and the mag- 
netic vector potential A, 



D hip _ 

= -V • u, 

Bt 

= -ciVlnp- 

Dt s ' 



(1) 



^ + V« + iVV- t( ) + / I (2) 
P P 



— = u x B - ?7/i J, 
at 



(3) 



where D/Di = d/dt + u ■ V is the advective derivative, 
B — V x A the magnetic field, J = V x Bj po the current 
density, and / is the random forcing function specified 
in B2001. The induction equation (y) implies a specific 
gauge for A; in Sect. |we discuss how to calculate gauge- 
independent magnetic helicity and helicity fluxes. Instead 
of the dynamical viscosity p (— const) we will in the fol- 
lowing refer to v = p/po, where po is the mean density in 
the domain (po = const owing to mass conservation). We 
use nondimcnsional units where c s = k\ = p = Mo = 1< 
Here, c s = const is the sound speed, fci the smallest 
wavenumber of the domain (so its size is 2ir), and po is 
the vacuum permeability. 

We use sixth-order finite difference and a third-order 
time stepping scheme. The number of meshpoints is usu- 
ally I20 3 , except for a few low Reynolds number runs. 

All variables are assumed to be periodic in the x- and 
y-directions. This implies that the vertical magnetic flux 
vanishes, because B z = d x A y — d y A x , so the integral of 
B z over x and y vanishes. Therefore the horizontal aver- 
age of the vertical field also vanishes. In the z-dircction 
we assume either also periodicity (if there is a halo), or 
we adopt a stress-free boundary condition where the field 
is purely vertical (if there is no halo). The latter states 
that d z u x — d z u y — u z = together with B x = B y = 0, 
and solenoidality then implies that d z B z = on the ver- 
tical boundaries. In terms of A, this boundary condition 
is d z A. x = d z A y — A z = 0. In models with a halo the 
extent in the z-direction is Air and / is multiplied by a 
mask function 



p(z) 



exp[ 




-z 2 /(H 2 -z 2 )} 



for 
for 



W < H, 
\z\ > H, 



(4) 



where H = n is the semi-height of the disc. In the models 
with halo we adopt periodic boundary conditions which 
has some conceptual advantages in that all the mag- 
netic energy and helicity that crosses the boundary at 
z = ±H can be accounted for. We also consider models 
where, in addition to the halo, we have imposed gravi- 
tational acceleration, so / — > / + g, where g = — V<f> 



Table 1. Summary of the main properties of the runs 
with vertical field boundary condition and with conduct- 
ing halo. The models Buoy 1 and Buoy 2 are stratified 
and allow for magnetic buoyancy effects. The parameter 
q = (B 2 ) I (B 2 ) gives the fractional magnetic energy in the 
mean field relative to the total magnetic field. 
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0.01 


0.01 
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Vert 2 


0.005 


0.005 
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Vert 3 


0.002 


0.002 




0.197 
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0.39 


Vert 4 


0.002 


0.001 




0.191 


0.166 


0.24 


Halo 1 


0.002 


0.0005 


0.02 


0.116 


0.137 


0.44 


Halo 2 


0.002 


0.001 


0.02 


0.118 


0.156 


0.60 


Buoy 1 


0.002 


0.001 


0.02 


0.124 


0.200 


0.61 


Buoy 2 


0.002 


0.001 


0.001 


0.147 


0.242 


0.71 



and $ = — <?oCOs(z/2) is a periodic gravitational poten- 
tial in \z\ < 2ir. This produces a density enhancement 
near the midplane at z = 0. The models presented be- 
low have go — 0.5, so the maximum density contrast is 
A lnp = 2go/Cg = 1. These models allow for additional 
flux loss via magnetic buoyancy. We first discuss the mod- 
els without halo and uniform background density. 

3. Description of the various runs 

The parameters of the various runs are summarized in 
Table § All runs develop a large scale magnetic field after 
some time. This large-scale field points in the horizontal 
directions and varies along z. However, in some cases there 
can be intermediate stages during which the large scale 
field is oriented differently. The different orientations of 
the mean field are best described by calculating the three 
different mean fields, whose energies are denoted in the 
following by K x , K y , and K z . The corresponding mean 
fields vary only in the x, y, and z directions, respectively, 
and are averaged in the two perpendicular directions; see 
also Eqs (8)-(10) of B2001. The mean field that varies 
only in the z-direction, and which is averaged over x and 
y, is denoted in the following by an overbar. 

3.1. Runs with vertical field boundary condition 

In Fig. [l| we show the resulting field structure near the end 
of run Vert 4 at t — 1150. Unlike the case with periodic 
boundaries where the presence of a large scale field was 
readily visible in two-dimensional slices of individual field 
components (see Fig. 4 in B2001), averaging over one of 
the two horizontal coordinate directions is required before 
the mean field becomes visible. 

Averaging over both horizontal coordinate directions 
clearly yields the profiles of the two non-vanishing field 
components, B x and B y ; see Fig. 0. Note that one of the 
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x average 




Fig. 1. Images of a cross-section and an x-average of B y . 
Vert 4, t = 1150. Dark and light indicate shades indicate 
negative and positive values, respectively. Note the high 
noise level in the cross-section. 
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Fig. 2. Horizontal averages of B x and B y . Vert 4. 



two components (here B x ) is approximately antisymmet- 
ric about the midplane, whilst the other (here B y ) is ap- 
proximately symmetric about the midplane and of larger 
amplitude. There are also some indications of boundary 
layer behaviour which was seen in highly nonlinear a 2 - 
dynamos (Meinel & Brandenburg 1990), with the same 
boundary conditions as in the present simulations, i.e. 

B X = By = 0. 

In Fig. H we plot the evolution of the magnetic en- 
ergy within the domain. The magnetic diffusion time is 
(rjk 2 ) -1 w 500. For t Ss 200, magnetic energy is mainly 
in the small scales, for which the helicity constraint is ir- 
relevant, so the initial exponential growth can occur on 
a dynamical time scale. At t — 400, the magnetic energy 
reaches a peak. At that time most of the large scale field 
energy is actually in K y , which denotes the energy of the 
mean field denned by averaging in the two directions per- 
pendicular to y, so this field varies only in the y-direction. 
This field is different from the anticipated field K x that 
is obtained by averaging in the x and y-directions. This 
K z field also builds up and saturates around t = 700, but 
it reaches a saturation level that is significantly less than 
in the case of periodic domains (B2001). Furthermore, as 
the magnetic Reynolds number increases, the fractional 
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Fig. 3. Evolution of kinetic and magnetic energy together 
with the magnetic energies contained in three differently 
averaged mean fields. Vert 3. 
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Fig. 4. Mean-field magnetic energy relative to the total 
magnetic energy, q = (B 2 )/(B 2 ), for different values of 
the magnetic Reynolds number, R m = u rms L/n, where 
L = 2tt is the scale of the simulation domain. 



energy of the mean field relative to the total magnetic 
field decreases; see Fig. [|. 

3.2. Runs with a conducting halo 

We have carried out runs with low and high halo con- 
ductivities, with and without imposed gravity. All these 
runs are similar with respect to the overall field evolution 
in that they all produce a clearly visible large scale field 
that grows slowly (on a resistive time scale), reaching fi- 
nally super-equipartition field strengths. These runs are 
therefore similar to those of periodic domains, but rather 
different to those with vertical field boundary conditions. 

In Fig. U we show a vertical slice of the resulting field 
structure for a run with gravity (go = 0.5) and a halo 
conductivity that is equal to the conductivity in the tur- 
bulent zone (in \z\ < 71"). Note that, in contrast to Fig. || 



4 



Axel Brandenburg and Wolfgang Dobler: Large scale dynamos with helicity loss through boundaries 





Fig. 5. Slices of the two horizontal field components for 
Buoy 2 with conducting halo and vertical gravity. The 
boundaries between halo and disc plane are indicated by 
white lines, t — 1400. Dark and light shades indicate neg- 
ative and positive values, respectively. 
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Fig. 6. Horizontal averages of B x and B y for Buoy 2. 
The boundary between halo and disc plane is indicated 
by the vertical dash-dotted lines, t = 1400. Note that at 
this late time significant amounts of magnetic flux have 
been diffused into the halo. 



the large scale field can clearly be seen even without av- 
eraging. The horizontally averaged field components are 
shown in Fig. |^. 

The evolution of kinetic and magnetic energies, as well 
as the magnetic energy of the large scale field, are shown 
in Fig. ^for Buoy 2. In addition to a run starting with just 
a weak seed magnetic field we also show how the energies 
evolve when restarting from a snapshot of Buoy 1 at t = 
800. 

In the following we study the results obtained in view 
of the helicity constraint, which was used extensively in 



B2001 in order to understand the slow growth of the large 
scale field and its final saturation level. 



4. Helicity constraint and helicity fluxes 

The resistively limited growth of large scale fields gener- 
ated by helical turbulence is primarily a consequence of 
helicity conservation. In the present case, however, there 
is an additional surface term which results from the flux of 
helicity passing through those boundaries, so the equation 
of magnetic helicity conservation becomes 



dH 

~dT 



= -2t]HqC - Q, 



(5) 



where H and C are magnetic and current helicities, respec- 
tively, and Q is the surface- integrated helicity flux through 
the boundaries. 

Owing to homogeneity of the system in the two hori- 
zontal directions we can consider integral quantities nor- 
malized per unit surface, so the current helicity (per unit 
surface) is therefore defined as 



C 



1 



L X Ly 



J BdV 



J ■ Bdz, 



(G) 



where the overbar denotes x and y integration and di- 
vision by the corresponding surface area L x L y (= 47r 2 ), 
i.e. a horizontal average. The magnetic helicity would be 
J A ■ B dz, but in this form it is not gauge-invariant and 
would change if a constant or some gradient field were 
added to A, which would leave B unchanged]^]. The gauge 
invariant helicity of Berger & Field (1984) is given by 



H 



(A + A )-(B-B )dz, 



(7) 



We recall that in a periodic domain, J A ■ B dV is auto- 
matically gauge invariant, because then there are no surface 
terms and J V<f> ■ B dV — — J <f)V ■ B dV = for any cf>. 
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where Bq = V x Aq is a potential field that satisfies 
V 2 A = 0, V • A Q = 0, A 0z = in z x < z < z 2 , (8) 
and 

A = A + V_l x (ipz), Vi^ = -B z onz = zi,z 2 ; (9) 

here, V_l = {d x ,d y ,0) is the horizontal nabla operator. 
Like A, we assume A and ip to be periodic. This implies 
that the horizontal average of V_l x (ipz) vanishes. This 
is the reason why we had to include an explicit A term 
in Eq. (S) which is not needed by Berger & Field (1984). 
Since the gauge-invariant magnetic helicity of the mean 
field plays an important role in the present context, we 
give an explicit derivation of this component of the mag- 
netic helicity in Appendix |a|. 

For the following analysis it is convenient to split the 
helicity terms into contributions from mean and fluctuat- 
ing fields, H = H me!m + H &uct , with [see Eq. (|A~5|)] 



Vert 4 



A ■ Bdz + z ■ (Ai x A 2 ), 



(10) 
z\ and Z2, 

(11) 



where A\ and A 2 are the values of A at z 
respectively, and 

iJ fluc t = / (a + a ) • (6 - b ) dz. 

J zi 

The latter expression is just Eq. (|^), but rewritten for the 
fluctuating fields, a = A — A, and correspondingly for 
the other variables. [We recall that Eq. ( |l(j| ) can also be 
written in a form similar to Eq. (fill) .] The components of 
the helicity flux, Q = Qmcan + Qfluct, are given by (see 
Appendix 



-(E 1 +E 2 ) 



Bdz 



\Z2 
Izi' 



and (see Berger & Field 1984) 
Qfluct = 2z ■ (e X ao) 
where 

E = r][M)J —u x B 



(12) 



(13) 



(14) 



is the electric field, and E = E + e is its decomposition 
into mean and fluctuating components. In Appendix^ we 
give the values of magnetic energy, magnetic helicity, inte- 
grated helicity flux, as well as the current helicity, for the 
eigenfunction of the a 2 -dynamo with the boundary condi- 
tion B x = By = 0. Next we present in Fig. ||the evolution 
of magnetic energy, magnetic and current helicities, and 
integrated helicity flux, for model Vert 4. Note that, in the 
steady state, the helicity flux is balanced by the current 
helicity term, i.e. 
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(15) 



Fig. 8. Evolution of magnetic energy M, magnetic helic- 
ity H, current helicity 2rjC and magnetic helicity flux Q, 
for mean (solid lines) and fluctuating (dashed lines) field 
components for Vert 4. (fiQ = 1 factors are omitted in the 
legends.) 



The small scale current helicity term, 277/10 C nuc t, may 
therefore be regarded as the source of magnetic helic- 
ity. Somewhat surprisingly, however, most of the helic- 
ity flux that is carried away results from the mean field, 
whilst most of the current helicity comes from the fluc- 
tuations. Similarly, |£f mC an| 3> |-fffluct|- This dominance is 
partly explained by the factor (k{/ki) 2 = 25 by which the 
large scales would exceed the small scales in the magnetic 
helicity relative to the current helicity [see B2001 after 
Eq. (41)]. 

It is then not too surprising, although perhaps some- 
what disappointing, that almost all the magnetic helicity 
flux is in the large scales. We note, however, that the frac- 
tional contribution of small scale helicity flux does increase 
somewhat with increasing magnetic Reynolds number; see 
Table |^, where quantities relevant for the magnetic helic- 
ity equation are given for the different models with ver- 
tical field boundary condition. The fact that the Q and 
C terms in Eq. (jl^) do not balance exactly is partly ex- 



6 



Axel Brandenburg and Wolfgang Dobler: Large scale dynamos with helicity loss through boundaries 



Table 2. Comparison of relevant quantities in the mag- 
netic helicity equation for different runs with vertical field 
boundary conditions. 
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Vert 3 
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mesh 
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0.0012 
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0.0012 


Q mean 


0.0007 


0.0017 


0.0017 


0.0010 


Mfluct /Afmcan 


0.3 
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1.6 


3.2 


-Cfluct / ^mcan 


-0.03 


-0.09 


-0.12 


-0.23 


Cfluct / Cmean 


-4 
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-26 
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Qfluct / Q mean 


-0.000 
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-0.007 


eh 


-2.5 
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plained by the finite value of dH/dt even towards the end 
of the simulation. A potentially important feature is the 
negative burst of small-scale magnetic helicity around the 
time when the small scale field saturates (i s» 300 in Figs || 
and |9|). In Halo 1 this is also associated with a similar neg- 
ative burst in Qfluct > as well as a with a period of enhanced 
small scale and large scale Poynting flux, 



ffluct = z-(exb) 



and P„ 



respectively. (We recall that in the case of the vertical field 
boundary condition P = because B is vertical.) 

Toward the end of the simulation with halo a signif- 
icant mean field has been accumulated in the halo (see 
Fig. ||). This is a somewhat unrealistic feature that re- 
sults from the use of periodic boundary conditions at the 
top of the halo. In reality the halo would extend further 
and field would be lost into the exterior (i.e. interstellar or 
intcrgalactic) medium. An additional artifact of the halo 
being saturated with a large scale field is the fact that 
at late times the Poynting flux is occasionally pointing 
inwards (Fig. which is unlikely to occur in reality. 

In Table || we also give the helicities and corresponding 
fluxes of the mean field in terms of the following nondi- 
mensional quantities: 



M ' 

ivi mcan 
Cmcan/^1 



Qmc 



(17) 



(18) 



(19) 



rma J1I mcan 



where k\ = 2it/L is the smallest wavenumber in the do- 
main (and usually the wavenumber of the large scale field) ; 



Table 3. Comparison of relevant quantities in the mag- 
netic helicity equation for different runs with different im- 
plementations of a conducting halo. Buoy 1 has a poorly 
conducting halo, whereas Buoy 2 has a halo conductivity 
equal to the disc conductivity. In Halo 1 the value of X] was 
so small that the full saturation field strength has not yet 
been achieved by the end of the simulation. The asterisks 
mark those quantities that are likely to change if the run 
was continued further. In all cases we use 120 x 120 x 240 
mesh points. 
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in our case fci = 1. We also have hq = 1, but we keep this 
factor in some of the expressions for clarity. 

Table | gives the same quantities for the model with 
a conducting halo. Here the overall helicity flux is gener- 
ally somewhat smaller than in the case with vertical field 
boundary conditions. However, the fractional contribution 
of small scale helicity flux is larger. The magnetic helicity 
production, which is dominated by 2?7/ioCfl uc t, is still sim- 
ilar in the two cases. In the following section we discuss 
why the losses of field with large scale helicity significantly 
lower the saturation field strength. 

5. Interpretation 

5.1. The helicity constraint 

In B2001 it was possible to understand the magnetic field 
evolution as a consequence of the magnetic helicity evo- 
lution. The current data show that most of the magnetic 
helicity and its corresponding flux is in the mean field. 
Thus, helicity conservation is governed approximately by 



, , -L* mean 

at 



-2rni (C n 



Cfluct) — Q 



(20) 



Furthermore, in B2001 it was found that for strongly he- 
lical fields both current and magnetic helicities are pro- 
portional to the large scale magnetic energy. The same is 
true here: from Table || we see that the ratio en is around 
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Fig. 9. Evolution of magnetic energy M, magnetic he- 
licity H, current helicity 2r/C, magnetic helicity flux Q, 
and magnetic energy (or Poynting) flux P, for mean (solid 
lines) and fluctuating (dashed lines) field components for 
Halo 1. All quantities are evaluated between the bound- 
aries z\ = — 7r and Z2 = n. After t w 400-500 the field 
in the halo becomes noticeable and effects of the periodic 
boundaries in the ^-direction begin to affect the evolution 
of the fluxes. 



2, and that £c is around 0.5-1, We also note that the he- 
licity flux is a certain fraction of the magnetic energy, i.e. 
eq fts 0.2-0.4. Using the fact that all these ratios are ap- 
proximately constant for the different runs, we may write 



n . ^mcan 

at 



-27/cfffcjMmcan + 2r]kl 



Cfl u 



where we have defined an effective magnetic diffusion co- 
efficient via 



2rj oS = {2r\E C + E Q k 1 1 u Tms )/e H . 
Equation (Efi) has the solution 



(22) 



AL 



V C flu 



1 



for t>t s , (23) 



2nefsk((t-t e 

where t s = — Aln(i?i n i/i? C q) is the saturation time of the 
total (mainly small scale) field, A is the kinematic growth 
rate of the dynamo, and i?i n i is the initial field strength. 
Equation ( |23] ) is analogous to Eq. (45) of B2001, but since 
r] e ff is much larger than 77, the time scale for reaching full 
saturation is reduced from (rikf)^ 1 to (r^ffkf)^ 1 due to the 
presence of boundaries. However, this result is obtained at 
the price of having reduced the saturation field amplitude. 
Therefore, all this flux term really does is cutting off the 
growth of the field earlier than otherwise, so at any time 
the energy is less than in the case of a periodic domain. 
Indeed, expanding the exponential for short enough times 
we find that 



Mr, 



t s ) for t > i s , 



(24) 



which shows that the initial growth time is still inversely 
proportional to the microscopic magnetic diffusivity. 

There is not much scope for increasing the rate of large 
scale field production in Eq. (pi}) . The value of Cfl uct does 
not strongly depend on the details of the model (cf. Ta- 
bles H and ||) . The only way out seems to be in lowering 
the value of en, which corresponds to making the field less 
helical. This is what happened when adding the effects of 
shear (Brandenburg et al. 2001). 

If r]/r] e s becomes very small the estimates given in 
Eqs. (^3|) and ( pjj ) underestimate the value of M mGan , be- 
cause there is always some mean field that results from im- 
perfect cancellations of the fluctuating field at the energy 
carrying wavenumber which is here kf. The resulting en- 
ergy of the mean field would be of the order (fci/fcf) 3 Mfl uct . 
In the present case we are however far away from this lower 
limit. 

In the presence of a conducting halo the helicity flux 
is smaller than in the case of the vertical field boundary 
condition; see Table || and Fig. |9| where we also show, in 
addition to magnetic helicity flux and current helicity, the 
evolution of magnetic energy, magnetic helicity, and the 
Poynting flux. Note in particular that <5 me an now fluctu- 
ates about zero. Consequently, there is no longer an ap- 
proximate balance between 2r]Ca uc t and Q m0 an- Instead, 
the Qmoan term becomes more important and would, in 
the limit of a periodic domain (Q — > 0), balance Cfl uc t, as 
was found in the case of periodic domains; see B2001. 

5.2. Comparison with a mean-field model 



In the case of a periodic domain it was possible to repro- 

(21) 

duce the slow saturation of the magnetic energy with an 
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Fig. 10. Normalized magnetic energy in one-dimensional 
mean-field models with periodic and open boundary con- 
ditions. The models with open boundaries allow for helic- 
ity flux out of the domain and yield significantly smaller 
field amplitudes when the magnetic Reynolds number, or 
X/r/kf, is large. In the periodic case the growth rate is 
A = a — ?7t, whilst in the case with open boundaries 
A = \{o? — r^)/rfx. For the periodic boundary condition, 
squares correspond to ao = 2 and triangles to a = 5. 
The data are obtained numerically using a time-stepping 
method. 



a 2 -dynamo with simultaneous a and r^-quenching of the 
form 



a 



1 + a B B 2 /B 2 



l + r lB B 2 /B 2 



(25) 



where as = t]b was assumed, and both quantities were 
proportional to the magnetic Reynolds number, as was 
suggested by Vainshtein & Cattaneo (1992). The reason 
this worked was because of the presence of microscopic 
magnetic diffusion which was not quenched, and which led 
to saturation such that [see Eqs (55) and (56) of B2001] 



B 2 



A 



aB ~cyr ~ ~~n ~ 0.01i? m (for a periodic domain), (26) -j 



where Bg n is the final magnitude of the mean field, A = 
ak\ — r^k 2 is the kinematic growth rate of the dynamo, 
Vt = V + f)t is the total magnetic diffusivity, and R m = 
UrmsL/r] is the magnetic Reynolds number with respect to 
the scale of the domain. 

In the present case of open boundaries there is an ad- 
ditional loss term and so Eq. (|26| ) no longer holds. A nu- 
merical integration of the nonlinear a 2 -dynamo equations 
with boundary conditions B x = B y = (Meinel & Bran- 
denburg 1990) shows that the term asB^ n / B 2 q saturates 
at a value of around 3 and is only weakly dependent on 
the value of ao; see Fig. where r\ + f? t o = 1- In the 
simulations, on the other hand, the final large-scale field 

— 1/2 

strength decreases with R m like R m , even though the 
total (mainly small scale) magnetic field strength is always 



close to the equipartition field strength; see Fig. ^ and Ta- 
ble PP. From this we can conclude that a B increases only 

1/2 

like i?m • This is to be contrasted with the case of peri- 
odic domains where Ba n /B eq was independent of R m and 
therefore «b a R m . It is important to realize, however, 
that the scaling of as with R m is secondary to having 
strong large-scale dynamo action. What really matters is 
the saturation level and the time scale on which full sat- 
uration is reached. These properties are controlled pri- 
marily by the magnetic helicity equation which, in turn, 
is sensitive to boundary conditions, anisotropies, and the 
presence of shear, for example (e.g. Brandenburg et al. 
2001). 

6. Conclusions 

The results presented in this paper have shown that open 
boundary conditions, and boundaries different from peri- 
odic ones, can have profound effects on the large scale field 
generation in helical flows. Firstly, unlike the case of peri- 
odic domains, a strictly force- free large scale magnetic field 
is no longer possible. Secondly, the loss of magnetic helic- 
ity through the boundaries severely limits the large scale 
field amplitude attainable in simulations. This is quite dif- 
ferent from the case of periodic domains, where the large 
scale field energy can exceed the kinetic energy, provided 
there is sufficient scale separation between the energy car- 
rying scale and the scale of the system (see Sects 3.5 and 
3.6 of B2001). For an open simulation volume with a verti- 
cal field boundary condition, however, the large scale field 
amplitude can be severely limited in a magnetic Reynolds 
number (R m ) dependent fashion. Again, this result can be 
modelled with a mean-field a 2 -dynamo, where the a-effect 
and turbulent magnetic diffusivity arc quenched. However, 
the i? m -dependence of these quenching expressions is here 
less extreme than in the case of periodic domains where 
the field can be force- free (B2001). This suggests that a 
dependence of a and rjt on B 2 alone is too simplistic, and 
that there might be an additional dependence on J x B 
or J ■ _B, for example, that would react specifically on non 
force-free magnetic fields. 

Our result that the quenching of a and r]t can be sen- 
sitive to external factors is not entirely unfamiliar. For 
example in the presence of large scale shear (correspond- 
ing to differential rotation) the quenching of a and 7y t was 
already found to be significantly reduced; see Branden- 
burg et al. (2001). A significant reduction of the quenching 
may explain why the cycle period in oscillatory dynamos 
does not need to scale linearly with the resistive time scale 
(which would be too long). Nevertheless, the growth time 
of a large-scale dynamo with helicity does seem to scale 
linearly with the resistive time scale in all cases investi- 
gated so far. 

In the simulations presented here the effects of open 
boundaries and magnetic buoyancy were too weak to pro- 
duce any sizeable small scale magnetic helicity fluxes. 
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However, this leaves unanswered the question whether 
small scale magnetic helicity fluxes are in principle ca- 
pable of reducing the long time scales that result from 
magnetic helicity conservation. In order to clarify this one 
may need to resort to models where small scale flux losses 
are somehow artificially enhanced. 

Finally, it should be mentioned that we have here ne- 
glected the effects of an equator where the sign of the he- 
licity would change. In the present framework this can be 
modelled by modulating the sign of the helicity of the forc- 
ing function in the ^-direction. Preliminary results have 
been reported in Brandenburg (2001b), but the results are 
quite similar to those in the case of open boundary con- 
ditions. Indeed an estimate similar to Eq. (0) has been 
obtained for the resulting magnetic energy of the mean 
field. 

It is now important to move on to more realistic global 
simulations which would automatically allow for helic- 
ity exchange across the equator and produce flux losses 
through winds and eruptions into the corona. It is indeed 
quite possible that significant helicity fluxes are only pro- 
duced by the combined action of stratification and shear 
(Berger & Ruzmaikin 2000, Vishniac & Cho 2000), which 
again would automatically be addressed by more realis- 
tic global simulations. Even if the magnetic helicity flux 
continues to be dominated by large scale contributions, 
it could be argued that such an effect could alleviate the 
helicity constraint, because saturation (albeit at a lower 
value) does occur earlier. It would seem necessary, how- 
ever, to compensate those additional losses by more pow- 
erful dynamo action. Again, the effects of shear may prove 
important, because it provides additional toroidal field 
that is non-helical and hence not subject to any constraint. 
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Appendix A: Gauge-independent magnetic 
helicity of the mean field 

In a periodic domain the magnetic helicity is gauge invari- 
ant, so AB (in z\ < z < z 2 ) can be made gauge invariant 
by linearly extrapolating A from z 2 to 2:3, where Aj, = A\ 
is assumed, and Ai are the values of A at z — Zi, where 
i = 1, 2, or 3. Linear extrapolation is used because this cor- 
responds to the simplest possible, current-free field. Thus, 
we can write for the magnetic helicity of the mean field 



A - Bdz + 



mean — / ' 5 - ' ' 1 - / * y - ±j y±Z, 

Zl J Z2 



A - Bdz 



where 
A 



A 3 - A 2 , 
A 2 + -i (z - z 2 ) 



Z3 - Z 2 

is the linear extension of A from z 2 to 23, and 
A 3 -A 2 



B=VxA=zx 



Z3 - z 2 



A B 



A3 - A 2 



z x 



23 - z 2 
Since A3 = A\ we have 



A 3 x A 2 

Z3 - z 2 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
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A ■ Bdz — z ■ (Ai x A 2 ), 



(A.5) 



which is independent of the choice of the value of z%. In 
particular, we may choose 23 = zi, in which case this 
definition is identical to Eq. (0) with Eq. (^) and Aq = A. 



Appendix B: Helicity flux of the mean field 

The evolution of the gauge-invariant magnetic helicity of 
the mean field is given by 



dt 



(A- B + A B)dz 



-z ■ (A 1 x A 2 + A x x A 2 ), 



(B.l) 
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where dots denote partial time derivatives, but, because 
A = — E in our gauge [cf. Eq. (0)], we have 

^#mea„ = -2^ " E ■ B dz - Z • (E X A) | J 

-2'(£ix3 2 +IixB 2 ). (B.2) 
The last two terms reduce to 

-z- [{Ex + E 2 ) x (A 2 - A x )] , (B.3) 
which can also be written as (i*^ + E2) ■ B dz. 

Appendix C: Helicity flux for the eigenfunction of 
the associated a 2 -dynamo 

The kinematic a 2 -dynamo equation with B = B(z,t) = 
{B x , By, 0) can be written in the form 

d t B = iad z B + r) T d*B, (C.l) 

where we have adopted complex notation, B = B x + \B y , 
and a and t/t are assumed constant. The eigenfunction 
that satisfies B = at z = ±tt is (e.g., Meinel & Branden- 
burg 1990) 

z ( cy \z \ 
B(z, t) = Bo cos - exp I \t , (C.2) 

2 v rn 27 

where Bq is a complex constant and A = \{p? — f^/rfr 
is the growth rate. We are interested in the marginally 
excited mode, A = 0, and so we assume for simplicity 
So = a = Tfr = 1, so we have 

B(z,t) = e- iz/2 cos(z/2). (C.3) 

The complex current density, J = id z B, is given by 

J(z,t) = ±exp(-iz), (C.4) 

and the complex vector potential is, apart from some ar- 
bitrary constant, 

A(z, t) — \ exp(— iz) — iiz, (C5) 

which satisfies d 2 z A = — J and B = id z A. For this field the 
magnetic energy M mean , the magnetic and current helici- 
ties £f mcan and C mea ii, as well as the integrated magnetic 
helicity flux Q mcan , are given by 

The Lorentz force is J x B = (0, 0, \ sin z), and the current 
helicity density is J ■ B = \ cos 2 (z/2), so 

so the field is neither fully force-free nor fully helical, but 
something in between. (For the mean fields of B2001 those 
values were and 1, respectively.) 



